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A formula to calculate a norm overlap between Hartree-Fock-Bogoliubov (HFB) 
states with the odd number parity (one quasi-particle excited states) is derived with 
, help of the Grassmann numbers and the Fermion coherent states. The final form 

(N 

of the formula is expressed in terms of a product of the Pfaffian for a neighboring 

o : 

■ even-even system (the zero quasi-particle state), and an extra factor consisting of 



(N 



the Bogoliubov transformation matrix and the anti-symmetric matrix in Thouless' 
HFB ansatz for the even-even system. 



INTRODUCTION 



^ ! The Hartree-Fock-Bogoliubov (HFB) method has been a powerful method in descriptions 

^ ' of nuclear states The reason is that the method can deal with the two most important 
^ correlations in interacting many-body nuclear systems, that is, deformation and pairing. 

These correlations are effectively taken into account by breaking relevant symmetries (the 
rotational and gauge symmetries). As a consequence, however, the conservation laws of 
angular momentum and the particle numbers are violated. 
' Restorations of these broken symmetries are achieved through quantum number projec- 

tions, but there have been difficulties to overcome in carrying out practical calculations of the 
projection. The major problem lies in the calculation of the so-called norm overlap kernels, 
which is necessary in the projections. The evaluation of the norm overlap is particularly dif- 
ficult in angular momentum projection because the rotational symmetry is associated with 
the non-Abelian SO (3) group. 
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An analytical formula was demonstrated to express the norm overlap by Onishi and 
Yoshida |2|. But due to a square root operation appearing in the formula, the relative sign 
(or phase) of HFB states at various points in the Euler space needs to be determined with 
respect to the reference HFB state, with an additional effort. 

Hara, Hayashi and Ring [3j were the first to perform a numerical calculation of angular 
momentum projection. They made use of the continuity and differentiability of the overlap, 
in order to determine the sign. 

However, it was later found that the assignment was sometimes very hard to be achieved 
due to a peculiar nature of the overlap. Such a case was seen in the cranked HFB wave 
functions, and it was discovered that the so<;alled "nodal lines" (a collection of zeros of 
the overlap) are 
presented in Ref. 



;he source of the problem |^. A method to overcome this problem was 
J], and an improvement to the method was recently found by the authors 



5|. With this method based on the Onishi formula, the sign problem was solved. 



Robledo proposed a totally different approach to the sign problem js]. Making use of the 
Grassmann numbers and the Fermion coherent state, he was successful to remove the square 
root operation in the norm overlap formula. His new formula is described by means of the 
Pfaffian. His approach to rely on the Grassmann algebra is not only mathematically elegant, 
but also quite powerful in practical computations of overlaps of many-body operators. There 
can be many applications to be discovered through his new methodology. 

In this paper, we would like to present such an application: a formula to evaluate a norm 
overlap between two HFB states with the odd number parity, which corresponds to nuclei 
with the odd-mass number. 

II. HARTREE-FOCK-BOGOLIUBOV STATES 

When the total number of constituent particles is even, the corresponding HFB ansatz 
(Thouless ansatz) is given as 

|HFB) = ATexp (j^ ) l^)- (1) 

Although this ansatz breaks the particle number conservation, the number parity is kept 
to be positive, or ( — 1)^" (n integer) In the above expression, the dimension of the 
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configuration space is given as M, and the (true) vacuum state |0) is defined as 

q|0)=0, (2) 
where the single-particle annihilation (creation) operator of the state i is expressed as q 

(4). 

Quasi-particle bases are introduced through a canonical transformation called the Bo- 
goliubov transformation W, which is 

' /3 





(3) 



where 



U V* 

W=\ |. (4) 

V u* 

Here, (3} and /3j are creation and annihilation operators for quasi-particles, respectively. The 
U and V are both M x M matrices, and can be regarded as the variational parameters 
in the HFB theory. These matrices satisfy the properties coming from the unitarity of the 
Bogoliubov transformation, that is, W^W = WW^ = X (X is the identity matrix) 

The anti-symmetric matrix Z appearing in the HFB ansatz is related to the UV matrices 

as 

z = {vu-^y. (5) 

The anti-symmetry of the above matrix can be easily confirmed by the properties possessed 
by the UV matrices. 

The normalization constant Af is calculated to be 



TV = VdetU, (6) 

by using the Onishi formula In many calculations of physical interest, such as one- 
dimensional cranked HFB states, U is a real matrix, so that A/" is a real number. 

The HFB state |$) is the vacuum of the quasi-particles, that is, the following condition 
is satisfied, 

m) = o. (7) 

In the quasi-particle bases, the corresponding energy spectrum is given as 

HnFB = Eo + Y,Ekl3ll3k. (8) 
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Excited states with many quasi-particles can be produced by operating the quasi-particle 
creation operators to the HFB vacuum. For example, a one-quasi-particle excited state is 
expressed in the framework of the HFB theory as Jj , 

\<^k) = Plm. (9) 

n 

The excited energy corresponding to this state is Ek- As explained in p. 250 of Ref. pLj], the 
new state has the negative number parity, which means that the state corresponds 
to an odd- mass nucleus, that is, a neighbor to the even-even nucleus described by |$). 

For the sake of simplicity, the isospin degree of freedom is not considered in the present 
paper, but the extension of the theory can be easily made. 

III. PFAFFIAN FORMULA TO NORM OVERLAP 

A. Case of even number parity 

The case of the even number parity was well studied by Robledo Q], and the Pfaffian 
formula of norm overlap kernels was derived for even-even nuclear systems for the first time 
in his work. But his conventions for the mathematical objects are slightly different from 
Ref. 8j, which is employed in the present work. For the sake of consistency in the following 
discussions, it will be convenient to derive the formula again here with our mathematical 
conventions. 

The essential point to derive the Pfaffian formula is to introduce the Fermion coherent 
state and its completeness. The Fermion coherent state 8| reads 

1^) =e-^'«''=^|0), (10) 

and it satisfies by definition the eigenvalue equation, 

cm = m, (11) 

where represents the Grassmann number. The Grassmann numbers follow the anticom- 
mutation rule, that is, ^i^j + ^j^i = 0. In the special case of i = j, there holds = 0. 
The completeness is given as 

I n^^*^^^^^p (-E^K^) = 1- (12) 
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Let us write a HFB state as the following: 

=f(Z(^'),ct)|0), (13) 

where an operator T is introduced as 

f(Z,ct)=exp|^iX]Z,,clc]j, (14) 

for p = 0,1. In this notation, T{Z, c^Y = T{—Z*,c), because Z is an anti-symmetric matrix. 
The overlap between two HFB states is thus expressed in the following way. 

($(0) !$(!)) (15) 
= (0|f(-Z(°)*,c)f(Z«,ct)|0) 

n«"e-^/'«^'^''(0|f(-ZW*,c)|O(^|T(Z«,ct)|0) 



In the last line, the operator T is replaced with a Grassmann-number quantity T due to 
Eq.( lTT|) . which is 



2 ^ ^ijCiCj J 
ij J 



T(z,r) = exp(^->^z,,e;e;i. (le) 

In addition, a property of 

(0|O = 1 (17) 

is used in the last line. 

Let us write the integrand in Eq. ffTSl) as G{^), that is, 

G{C) = e-^/3«^^^T(-Z(°)*,0T(zW,r) (18) 

As demonstrated by Robledo [6] , G can be summarized to be a Gaussian with the Grassmann 
number: 

G(0 = exp (Ic'^c) , (19) 



2 

where Grassmann vectors in his work are defined as 



C* = (^1;'C25 ■ ■ ■ 5^/56; ^2, ■ ■ ■ ,^m), (20) 



and a matrix Z with the 2M dimension is equal to 

Z = . (21) 

X is the M X M identity matrix. Apparently, Z is anti-symmetric. 

In the present work, the following ordering for Grassmann vectors is employed, 

= i^h^h ■ ■ ■ ,Cm,^m,^m-i, ■ ■ ■ ,6), (22) 



which is the same as in Ref. 



a 



The transformation from C to is achieved as 

^=LC, (23) 



where a linear transformation is given by 

L=l . (24) 

yo A 

A M X M matrix A is defined as 

Ajj = 5i+j^M+i- (25) 
This matrix satisfies L = = L^^ (A = A* = A^^), and 

det(L) = det(A) = ^-1)M(m-i)/2_ ^26) 

After the transformation, the Gaussian G is rewritten as 

G(C) = G(0 = exp . (27) 

The relation between Z and X is X = L*ZL, so that 

/ -A \ 

X = . (28) 

y A -AZW*A j 

Thanks to a mathematical theorem j^, [lo|, it is always possible to find a decomposition 
of an anti-symmetric matrix into a matrix product 

X = i?*Ji?, (29) 

where i? is a regular matrix (3i?^^) and J corresponds to a canonical form of X, that is. 



O X 
-X O 



(30) 
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The new Grassmann bases t] associated with the canonical form is obtained by a hnear 
transformation of the original bases ^, 

V = (31) 

For the sake of convenience in subsequent discussions, let us write the inverse transformation 
as the following 

^.n-^,^h-]J^^^^A(A. (32) 

It should be noted that the index of rj and ^ runs in a reverse order, as in Eg. (1221) . 
In the new bases fj, G has an expression of 

G(0 ^ G{v) = exp (^r^\R-yXR-'ri\ (33) 



1 / 

exp(-r/*Jr7) = exp f ^r/* 



f]M+l-a 



Let us introduce a notation a here for convenience in subsequent discussions, which is defined 

as 

a = M + l-a. (34) 

Using the property of the Grassmann number [rif = 0), a Taylor expansion of the exponential 
can be greatly simplified to a sum of bilinear polynomials. It is thus possible to write G as 

M 

G{n)=l[{l + v:v^). (35) 

a 

Furthermore, it is important to understand that only the non- vanishing contribution of such 
polynomials in the 2M-dimensional Grassmann integral comes from the integrand in which 
all the bilinear pairs {f]ayVa) appear. After an expansion of the above product, such paired 
term appear in the form of Yla^aVa- To convert this expression to the standard paired 
expression, there is a useful identity 

M M 
a a 

where M is an even integer. 

With this identity, the integral of the paired polynomial is simplified, and the calculation 
is easily done. 

/M M 
U^VldVall^lv^ = ■ (37) 

a j3 
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The integral of the other polynomials give a null contribution to the overlap calculation of 
our interest. 

Noting that the Jacobian for the bases transformation R is det(i?), the integral of G{^), 
or the overlap, becomes 

($(o)|$(i)) = / \[di^di^G{i) (38) 
= (-1)^^/Met(i?) f l[dv:dVo.l['^;v^ 

''a 13 

= (-l)^/2(_l)^det(i?). 
Using another identities, that is, 

Pf(X) = det(i?)Pf(JJ), (39) 

and 

Pf(J) = (_l)M(M-l)/2^ 

the final expression is obtained as 

^$(0)|^{l)^ _ (_^)M(Af+2)/2pf^X) _ pf^X)_ (4^) 

The phase factor gives rise to (— 1)*^(*^+2)/2 = _|_]^ for even M, and it is different from 
Robledo's formula, which is (— This is because of the difference in the definitions 
of the Grassmann vectors, Eqs. (12(1 and (12^ . and it is simply explained as Pf(X) = 
det(L)Pf(Z) = (-l)AnA/+i)/2pf^^) i_ 



B. Case of odd number parity 

Next, let us consider two HFB states with the odd number parity. Following Eq.([9]), they 
are expressed as 

|$(o)^^^t(o)|^(o)^. !$(!))= /3t«|$(i)). (42) 

Here, the HFB states |<I>(p)) {p = 0,1) are considered to be states with the even number 
parity that are given by Eq. (IT^ . 



1 (_i)M(M+i)/2 ^ (_i)Af(A/-i)/2 because (-1)^ = 1 for M being even. 
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Below, a formula is derived for an overlap between the two states with the odd number 
parity, that is, 

{^'^W) = {<^''Wpr\<^'''), (43) 

which is expressed in terms of the Bogoliubov transformation matrices {U and V) and the 
inverse of X given in Eq. (1281) . 



It may be worth noting that a re- arrangement of a product = (3^^^ + 

(^t(o)f/(i) _j_yt(o)y(i)^^^, does not help very much to simplify the overlap because the vacuum 
condition is applied only to the associated annihilation operator. In other words, /Sj^^l^'^'^^) = 
only when p = q. 

In the single-particle bases, the product Pk'^ PI^^^ is expanded as 

/3f = 6,,,(c, ct) + kkk'ic, ct). (44) 

Two operators & and A are bilinear functions of the single-particle creation and annihilation 
operators, which are defined as 

6^^' = E ('^f ^^S' + '^''4cj) , (45) 

^^^y = E {^'f^^^^ + ^f^4) ■ (46) 

Let us call & the normal operator while h the dangerous operator. The matrix elements in 
the right hand side are given in terms of the UV matrices. 

, ( '^r s^f\^( (vt)i«)t7<;> ([/.)£)[/<;) \ 
cry W)^yl^ mf^vff) 

The calculation of the overlap given in Eq. (l43l) is thus reduced to a sum of the overlaps of 
the normal and dangerous operators, with respect to the HFB states with the even number 
parity, that is, 

= {^^'^\&kk'\^^'^) + (48) 
In the following subsections, each term in the right-hand side is considered separately. 

C. The normal operator & 

The essential ingredient of the calculation here is an evaluation of ($*^°^|qc|,|$*-^''). We 
begin with an insertion of the completeness for the Fermion coherent state between the 
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creation and annihilation operators as 

(<I>(°)|qc]|<I>(i)) (49) 

, M 



M 



a 

„ M 

^ J U'^Cd^-^mCr (50) 



The G{^) has the common structure seen in the even number parity case, that is, Eq.f lT9|l . 

Considering that ^j^* gives rise to three kinds of terms consisting of bihnear expressions 
of rja and r^*, what we need to calculate in Eq.f lSOj) are G{ff)rjirjj, G{f])ri*ri*, and G{ff)rj*rjj. 
However, from a simple analysis, the first two cases go to zero after integrations. It is thus 
enough to consider the last case. 

Noting the reverse order in the index for rj and ^, we have 

M 

C = + (51) 

i=i 

M 

i^ = E(^2l),,^; + (7^22),,^J. (52) 



From the property of the Grassmann integral, which is given in Eq. fl37j) . the non- vanishing 
contribution comes from the term of Tj-jT]*,. Due to an identity relation shown in Eq.( l37|) . 
there must hold a relation between the indices j and j', which is to be explained below. 

Let us consider a product (1 + ri*,rij,){l + rij,rij'){l + ri*rij){l + 'rj-.rjj)rj-jrj*i, assuming 
that j 7^ j' . The first four factors are always included in a representation of G{f]), 
given in Eq. fl35|) . Not only the factors commute mutually, but also with the other fac- 
tors in Eq. (l35|) . Then, a re-ordering of the product simplifies the first and third factors, 
thanks to a property of the Grassmann number, that is, (1 + rj*,rj-ji)rj*,(l + rj%,rjji){l + 
VjVj)Vj{^ + rjjrij) = + ?7j,?7j/)?7j(l + rjjrjj) Then, an expansion of the product gives 

rise to ri*,rij — rj*,rjjrj-.rj-j — ri*,rij/rij,rij + rj*ir]jirj-.rj-jrj%,rjj. All these four terms give a null con- 
tribution to the integral due to a presence of unpaired Grassmann numbers. In order to 
maintain the pair structure as seen in Eq. (l37|) . there must hold j = j'. Then, the relevant 
product becomes (1 + ri*rij){l + rijrij)rj*rjj = ri*rij{l + rj-.rjj) = —Tj--rj-jrj*rjj, which does not 
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vanish after integration. The result is thus summarized in an identity 

„ M 

/ n dvldiiMvHVk^ = (53) 

Noting that the Jacobian is given by det(_R) and putting all the above result together, 
the whole integral Eq. flSOj) becomes 

HdCd^^GmCj (54) 

a 

det{R) [ HdT^ldvMn) J2 fe^7 (^2i),^ (7ei2) 



J7 

/37 



-V*V^ (^22)i^ (^ll)j7, 

= (-i)3*^/Met(i?) (7^2l7^i2 -7^227^*0,^. 
= Pf(x) (7^2l7^i2-7^227^i0,, 

In the last line, we have used relations Pf(X) = det(i?) Pf(J) and Pf(J) = (-1)^^(a^-i)/2. 
Also, a fact that M{M + 2)/2 is an even integer was used. 

The similar calculation is performed for ($'^°^|c|cj|$''^^). In this case, it is important 
to exchange the order of a product of the creation and annihilation operators, that is, 
Cjc] + 6ij, for the convenience in applying the completeness of the Fermion coherent 



c c 



state. The first term corresponds to the exactly same result as obtained above, except the 
sign and the change in the indices {i ^ j). Whereas, the second term is a C-number, so 
that the overlap is simply proportional to ($^°) |$*^^^). 

The final expression for the normal operator becomes the following 

($(°|6fcfc'|$(')) = ($(°)|$(^))Tr (t'^'^'^i + S'''''^2) , (55) 

where 

mi = 7^l27^*l - tZuTzI^, (56) 

m2 = 7^227^*n - 7^2l7^*2 + x. (57) 

In obtaining the above expression, the result obtained in Eq. lHTi) is also used. 

D. The dangerous operators R 

The essential ingredients in this subsection is overlaps of the so-called "dangerous terms" 
in the HFB theory, which are expressed as ($''*^^|cjcj|$*^^'') and {^^^''\ciCj\^^^'>) . Because a 
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complex conjugate of one type of the dangerous terms corresponds to the other, it is sufficient 
to consider a mathematical analysis for one of the two terms. Let us take the first type here. 
That is, 

($(°)|4c;i$W) (58) 



M 



n dCd^^e- ($(0) I ,t,t I $(1) ) 

a 
M 

a 

„ M 

= U^maGmCy (59) 

According to Eqs. flFI]) and flS2]) . a product ^*^* is expanded in terms of bilinear polynomials 
of 1] and 1]*. As discussed in the previous section, only the type of terms rjlrn^, contributes 
to the integral if k = k'. The result is given in Eg. (1531) . The final result is thus obtained as 

($(o)|etct|$(i)) = Pf(X) (7^n7^l2 - 7^l27^10,, ■ (60) 

The other dangerous term can be obtained in a similar way, 

($(°)|qc,|$(i)) = Pf(x) (7^2l7^*2 - 7^227^*2l),,■ • (ei) 

By combining these results, the dangerous part becomes 

= (<l>(°)|<l>(i))Tr (^IC'^'^'Di + r'^'^'Da) , (62) 



where 



Di = 7^227^2l - 7^2l7^*22, (63) 

2)2 = 7^l27^1l - 7^n7^*2. (64) 



E. The final expression 

The sum of the overlaps of the normal and dangerous operators gives rise to the final 
form of the overlap formula in the case of the odd number parity. 

Before writing down the formula, however, it is worth considering one more thing here. 
In the above discussion, we introduced the quantities expressed in terms of the inverse of 
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the transformation matrix i?, that is, OTj and Di {i = 1,2). It is possible to express these 
quantities directly through the inverse of X. The demonstration can be shown simply by 
taking the inverse of the both sides of Eq.f l29|) . The result is 

= -R~^I{R~y (65) 

'^'V (66) 

The formula is therefore obtained as 

($f = ($W|$«)Tr(W^'='X-^ + 5'=^'). (67) 

The final form of the formula is independent of the bases transformation matrix R. 

The advantage of this formula is that the overlap of the odd-number-parity is expressed 
in terms of the quantities obtained for the neighboring even-even nucleus l^*^*)): the overlap 
of the even-number-parity, i.e., Pf(X), and the Bogoliubov transformation W. Although 
the inverse of the matrix X needs to be computed for the formula, the matrix itself can 
be expressed through the quantities calculated for the even-even system (See.Eq( 12S]) ). In 
other words, the quantum number projections can be done simultaneously for an even-even 
nucleus and the neighboring odd system, without significant efforts. 

The similar procedure can be applied to a derivation of formulae for multiple quasi- 
particle excited states (with more than one quasi-particle). The basic structure is a product 
of the Pfaffian of a certain even-even system and the factors representing each quasi-particles 
expressed in terms of W, including the inverse of X. 



IV. SUMMARY 



A formula Eq. fl67p was demonstrated so as to calculate a norm overlap for the HFB states 
with the negative number parity (one quasi-particle states), which correspond to odd- mass 
nuclei. The Grassmann algebra and the Fermion coherent state are employed, so as to allow 
the Pfaffian to describe the overlap. 

The formula has a factorized structure, which consists of the norm overlap for an even- 
even system and part described in terms of the Bogoliubov transformation matrix, as well 
as the inverse of the matrix given in Eq.f l2S]) . This structure is beneficial because both of 
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the systems with the positive and negative number parities can be studied at the same time 
by means of angular momentum projection. 

Recently, Avez and Bender presented the similar work 



11| . as well as Bertsch and Robledo 



12| . These works, including ours, result in the Pfaffian, which were initially demonstrated by 



Robledo. However, in their HFB wave functions, fully blocked unpaired particles are assumed 
{V = 1,U = in terms of the Bogoliubov transformation), which are different from our 
ansatz Eq.([n]) based on the UV exchange approximation for multiple quasi-particle excited 
states. In addition, the mathematical representation of the final results are significantly 
different from each other. 



Acknowledgments 

The authors thank for useful discussions with Professor R. C. Johnson. This work is 
financially supported with a research grant in Senshu University. 



[1] P. Ring and P. Schuck, Nuclear Many-Body Problem, Springer- Verlarg, 1980. 

[2] N. Onishi and S. Yoshida, Nucl. Phys. 80, 367 (1966). 

[3] K. Kara, Y. Hayashi, P. Ring, Nucl. Phys. A358, 14 (1982). 

[4] M. Oi and N. Tajima, Phys. Lett. B 606, 43 (2005). 

[5] M. Oi and T. Mizusaki, in preparation. 

[6] L. M. Robledo, Phys. Rev. C. 79, 021302R (2009). 

[7] N. Onishi and T. Horibata, Prog. Theor. Phys. 64, 1650 (1980). 

[8] J. W. Negele and H. Orland, Quantum Theory of Finite Systems (MIT Press, Cambridge, 

MA/London, 1985). 

[9] B. Zumino, J. Math. Phys. 3, 1055 (1962). 

[10] P. Benner, R. Byers, H. Fassbender, V. Mehrmann, D. Watkins, Electronic Transactions on 

Numerical Analysis 11, 85 (2000). 

[11] B. Avez and M. Bender, arXiv[nucl-th]:1109.2078vl (2011). 

[12] G.F. Bertsch and L. M. Robledo, arXiv[nucl-th]:1108.5479vl (2011). 



